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ABSTRACT 

The first structures were proto-voids formed in the primordial plasma. Vis- 
cous and weak turbulence forces balanced gravitational forces when the scale 
of causal connection Lh = ct ^ Lsv = {l^/pGY^'^ ~ Lst = e^^'^/{pG)^^'^ at 
time t ~ 30,000 years (10^^ s), with c the speed of light, 7 ~ 1/t the rate of 
strain, and u the kinematic viscosity, where Lsv and Lst are viscous and tur- 
bulent Schwarz scales of hydro-gravitational theory (Gibson 1996). The photon 
viscosity u ^ A x 10^^ m^ s~^ allows only weak turbulence from the Reynolds 
number Ren = c^t/v ^ 200, with fragmentation at pL^gy IO^^Mq to give 
proto-supercluster voids, buoyancy forces, fossil vorticity turbulence, and strong 
sonic damping. The expanding, cooling plasma continued fragmentation to proto- 
galaxy-mass ~ IO^^Mq, with p ^ lO"^"^ kg m~^ and 7 i=a 10~^^ s~^ preserved as 
fossils of the weak turbulence and first structure. Turbulence fossilization by self- 
gravitational buoyancy explains the 5T/T pa 10~^ cosmic microwave background 
temperature fluctuations, not sonic oscillations in cold-dark-matter fragments. 
After plasma to gas transition at t ~ 300,000 years (10^^ s), gas fragmenta- 
tion occurred within proto-galaxies at Lj ^ lO^Lsv and Lsv ~ Lst scales 
to form proto-globular-star-cluster (PGCs) clouds of 10^^ small-planetary-mass 
primordial-fog-particles (PFPs). Dark PGC clumps of frozen PFPs persist as 
inner-galaxy-halo dark matter, supporting Schild's 1996 quasar-microlensing in- 
terpretation. Non-baryonic dark matter, with D ^ 10^^ m^ s~^, diffused into the 
plasma proto-cluster- voids and later fragmented as outer-galaxy-halos at diffusive 
Schwarz scales Lsd = (-D^/pG)^/^, indicating m «i 10~^^ kg weakly-coUisional 
fluid particles. Observations (Gibson & Schild 2003) support the theory. 
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1. Introduction 

We consider the hydrodynamic evolution of the hot big bang expanding universe after 
mass-energy equality to determine when gravitational forces were first able to form struc- 
ture under the influence of viscous and turbulent forces. All flows of plasmas and gases 
with large Re = Sv x L/u > Recr are unstable to the formation of turbulence accord- 
ing to the 1883 Reynolds number criterion for transition, where Recr is a finite critical 
value of Re above which laminar flows are impossible, 6v is the velocity difference on scale 
L and u is the kinematic viscosity of the fluid. From the flrst universal similarity hy- 
pothesis for turbulence (Kolmogorov 1941), the universal critical Reynolds number value 
Recr ~ 25 — 100 apphes to the Hubble flow as it does for all others. Gluon-neutrino- 
photon viscosity values u ~ cH before mass-energy equality at 25,000 years give subcritical 
Re ^ 1. Cosmic microwave background anisotropy extended Self Similarity (ESS) coeffi- 
cients (Bershadskii & Sreenivasan 2002) closely matching those for high Reynolds number 
turbulence arc attributed to fossils of pre-inflationary turbulence (Gibson 2001) at tempera- 
tures 10^^ — 10'^^ K too hot for large lepton viscosities to exist that might otherwise prevent 
big bang turbulence. Predictions of spectral forms and other flrst and second order turbulent 
flow parameters from Kolmogorovian universal similarity theories for turbulence and turbu- 
lent mixing have been widely validated in numerous atmospheric, oceanic and laboratory 
flows and numerous fluids (Gibson 1991). No experimental counterexamples exist, either for 
the Reynolds number turbulence transition criterion or Kolmogorovian universal similarity 
at low order. Linear stability theories suggesting the possibility of steady inviscid flows 
(Rayleigh 1880) have been recognized as unreliable for real fluids since the (Prandtl 1921) 
discovery of viscous instabilities and because much larger values of Recr are predicted by such 
theories than observed in laboratory experiments and numerical simulations (White 1991). 

Prior to the 1993 discovery that the anisotropics 5T/T of the cosmic microwave back- 
ground temperature are very small (~ 10~^) it was consistently assumed by all authors that 
Re values of the expanding universe would be supercritical (^ 100), so that both the plasma 
and the subsequent gas would be strongly turbulent with primordial turbulence the crucial 
factor in all subsequent gravitational structure formation. Density fluctuations produced and 
mixed by the turbulence would trigger gravitational collapse to form structures such as stars 
and galaxies at mass scales determined by the primordial turbulence. Prom energy argu- 
ments vol Weizsacker 1951 showed the Jeans 1902 hnear acoustical criterion for gravitational 
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instability in ideal fluids fails in strongly turbulent flows. He proposed that Kolmogorov's 
incompressible turbulence expression 5v ~ L^/^ for velocity differences 5v between points 
separated by distances L should be used to compute the turbulent kinetic energy of a pos- 
sibly unstable gas or plasma cloud, asserting that the turbulent kinetic energy of the cloud 
should be less than the gravitational potential energy as the criterion for gravitational in- 
stability in such clouds. Chandrasekhar 1951 also rejected the Jeans 1902 criterion for the 
gravitational instability of strongly turbulent flows but overlooked Kolmogorov's theory in 
any form and simply added a turbulence pressure pt ~ p{5vY to the fluid pressure p in the 
expression for Jeans's length scale Lj 

Lj = Vs/{pGf"^{p/p'Gfl\ (1) 

where p is the density, G is Newton's gravitational constant and Vs is the sound speed, to 
give a Chandrasekhar turbulent Jeans scale Ljct = [(p + Pt)/p^G']^/^. Star formation rates 
in the cold molecular clouds of the Galaxy disk are about 50 times less than expected from 
Eq. 1, presumably because Ljct > L > Lj, where L is the size of the cloud, Schefiler and 
Elsasser 1988, p438. Doppler broadened molecular absorption lines give strong evidence of 
Kolmogorovian turbulence in such clouds (Falgarone and Phillips 1990). A dissipation rate 
e ~ 10~^ m^ s~^ is estimated from the third order velocity structure function measured in the 
Ursa Major cirrus cloud (Miville-Deschenes et al. 1999), giving Lst ~ 8 x 10^^ m (p pa 10~^^ 
kgm~^) much larger than the cloud size so that star formation is prevented (see Table 1). 

The hydro-gravitational theory (HGT) of gravitational structure formation (Gibson 1996) 
abandons the Jeans 1902 theory in its entirety; not only for strongly turbulent flows but for 
flows that are weakly turbulent or nonturbulent. Lj in Eq. 1 should not be interpreted as 
either the minimum scale or maximum scale of gravitational instability as proposed by Jeans 
1902. Such misinterpretations have resulted in the dark matter paradox. For a self gravitat- 
ing ideal gas of nearly uniform density p and temperature T, Lj represents the maximum 
scale of acoustical pressure and temperature equilibration L/c = {RT/ pGY^'^, where it! is 
the gas constant and p is the pressure (Gibson and Schild 1999ab). Such a field of nearly imi- 
form plasma with known properties formed after the big bang and turned to gas at 300,000 
years. From HGT, non-acoustic density perturbations in the primordial plasma and gas are 
absolutely unstable to structure formation, and viscous or weakly turbulent fluid forces at 
Lst ~ Lsv ~ Lk = {v^/ey^^ Lh = ct, or diffusion at Lsd determine the smallest scales 
of gravitational instability, not Lj, where L^ is the Kolmogorov scale and Lh is the Hubble 
scale of causal connection. In the hot primordial plasma Lj > Lh, so by the Jeans 1902 
criterion no structure could form. Cold- dark-matter (CDM) non-baryonic fluid was invented 
with small Lj values to permit gravitational structure formation consistent with observa- 
tions (Padmanabhan 1993). However, the necessarily strong diffusivity Dcdm ^ cH of the 



-4- 



weakly coUisional non-baryonic dark matter in the plasma epoch prevents its condensation 
and rules out CDM models (Gibson 2000) because {Lsd)cdm ^ Lh in the plasma epoch. 

To correct the Chandrasekhar 1951 expression, the turbulent pressure ~ p{5vY should 
be substituted rather than added to p in Eq. 1 and the complete Kolmogorov 1941 expression 
5v ^ {eLy^^ should be substituted for Sv. Solving for the critical length scale at which 
inertial forces match gravitational forces gives 

LsT ^ e'l'l{pGfl\ (2) 

where Lst is defined as the turbulent Schwarz scale (Gibson 1996) and e is the viscous 
dissipation rate of the turbulence. 

If the turbulence of the primordial plasma flow is weak, as indicated by the small CMB 
fluctuations, then viscous forces Fy pv^l? determine the smallest scale of gravitational 
instability, balancing gravitational forces Fq ~ p^GL^ at the viscous Schwarz scale Lgv-, 
where 

Lsv={vilpGfl\ (3) 

V is the kinematic viscosity of the fluid, 7 is the rate of strain, and p is the density. The 
turbulent Schwarz scale of Eq. 2 is closely related to the Ozmidov length scale = 
[ejl^^yl'^ of stably stratified turbulent flows, where the stratification frequency has a 
physical significance similar to the inverse free fall time (pG)^/^ and is derived by matching 
turbulence forces with buoyancy forces to find the critical length scale. The viscous Schwarz 
scale of Eq. 3 near Rtcr is analogous to the buoyancy-inertial- viscous scale Lbiv = (z^/iV)^/^ 
that arises in fossil turbulence theory (Gibson 1999a). Turbulence is strongly inhibited 
and rapidly fossilized by buoyancy forces in the ocean and atmosphere at Lr scales, and 
astrophysical turbulence is strongly inhibited and fossilized at Lst scales in self gravitating 
fluids. Because kinetic and gravitational forces of a flat universe arc closely matched at 
the horizon scale Lh, it follows that whatever turbulence levels existed at the time of first 
structure formation (when, for the first time, Lsv ~ Lst < Lh) would be rapidly damped 
by buoyancy forces and the horizon length, density, mass, and the hydrodynamic parameter 
{e or 7) preserved by hydrodynamic fossils. 

Silk and Ames 1972 suggest that the large size of Lj ^ Lh = ct in the plasma epoch with 
sound speed Vs ~ c/S^^"^ prevents gravitational condensation of plasma by the Jeans 1902 
criterion. By their galaxy formation theory, strong turbulence produced density fiuctuations 
that served as nuclei for galaxy formation at the time of photon decoupling when the sound 
speed Vs dramatically decreased by a factor of 3 x 10^. Other studies claiming that strong 
primordial turbulence should set the scale of galaxies include Gamov 1952, Ozernoi and 
Chernin 1968, Ozernoi and Chernin 1969, Oort 1970, and Ozernoi and Chebyshev 1971. 
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All such strong turbulence theories of structure formation were rendered moot by the 
1993 measurements of very small temperature fluctuations 5T /T ^ 10~^ ~ 5v/v ^ 5p/ p ~ 
5p/p ~ 5a/ a in the cosmic microwave radiation (CMB) data from the 1989 COsmic Back- 
ground Explorer (COBE) satellite, rather than values oi5v/v ~ 10^^ — 10"^ that would result 
from fully developed turbulence, where a is the cosmic scale factor and 6{T,v, p,p,a) repre- 
sent fluctuation magnitudes. A subcritical horizon scale Reynolds number Ren = cH/u < 10 
at the time 10^^ s of plasma-gas transition requires an enormous kinematic viscosity u > 10'^^ 
m^ s"^ to be subcritical, much larger than u fa 10^^ m^ s""^ estimated for the primor- 
dial plasma then (Gibson 2000). For a terrestrial comparison, the kinematic viscosity of the 
Earth's upper mantle is u 10^^ m^ s~^ from glacial rebound rates (Professor Robert Parker 
of SIO, personal communication). Implicitly it has been assumed in the astrophysics liter- 
ature after these COBE observations that the Hubble flow of the expanding universe must 
somehow be intrinsically stable to turbulence formation, independent of Reynolds num- 
ber. Textbooks on structure formation in the universe such as Padmanabhan 1993 make no 
mention of viscosity, diffusivity, turbulence, or Reynolds number in their discussions of the 
process. No reference in the literature has been found that attempts to justify this implicit 
(and unwarrented) assumption. An example of strong turbulence generated by the Hubble 
flow is shown in Figure 1. Powerful Hubble flow drag forces separate protosuperclusters, 
protoclusters, and protogalaxies as they form by gravitational fragmentation in the primor- 
dial plasma and early gas epochs according to HGT. Hubble flow galaxy Reynolds numbers 
of order 10^^ shown in Figure 1 have decreased to values « 10^ or less at present. 

The assumption made by CDM hierarchical clustering cosmology models (CDMHCCs) 
that the Hubble flow is stable to the formation of turbulence is inconsistent with the universal 
similarity theory of turbulence, which is the basis of HGT. Strong turbulence in the plasma 
epoch is ruled out from HGT by the small values 6T/T ^ 10^^ of the CMB observations, so 
buoyancy forces resulting from gravitational structure formation must have dominated the 
damping of turbulence because viscous forces are inadequate and no other fluid forces exist. 
Turbulent transition cannot fail by lack of triggering perturbations since 5T/T fluctuations 
are observed at scales L > ctin the CMB that can nucleate growth of vorticity and structure 
once they enter the horizon. Neither can it be argued that a lack of time prevents self 
gravitational or fluid mechanical nonlinearity. Once Lx < Lh at turbulence transition 
the eddy overturn time is t. The decreasing viscous stresses in the baryonic component 
permit fragmentation of superclustcr to galaxy masses in the plasma epoch so that the 
hierarchical clustering of subgalactic scale CDM halos to form these structures in the gas 
epoch is unnecessary, even if such small CDM halos were physically possible (they are not). 
Because the nonbaryonic dark matter is necessarily strongly diffusive, such small CDM 
halos are excluded by HGT (Gibson 2000). Observations of galaxy-QSO correlations and 



-6- 



discordant cluster red shifts rule out CDMHCCs (Gibson & Schild 2003). CDMHCCs are 
also excluded by observed density distributions near galaxy cluster cores that fail to match 
universal forms computed by numerical simulation (Sand et al. 2002). 

In the following §2 we consider whether an inviscid expanding universe is stable or 
unstable to the formation of turbulence. If it is unstable according to the conventional 
Reynolds number criterion, what constraints on viscosities and structure formation in the 
plasma epoch can be inferred from observed CMB anisotropics? We then examine the 
hydrodynamic parameters and structures to be expected from the (Gibson 1996) nonlinear 
gravitational structure formation theory during the plasma epoch, in §3, and in the early 
gas epoch, in §4. Gonclusions are summarized in §5. 

2. The absolute instability of inviscid flows 

The instability of expanding flows is discussed in §23 of Landau and Lifshitz 1959. The 
equations of momentum conservation in a fluid may be written 

dv/dt^ -VB + vxuj + vV'^v + FM + ... (4) 

where B = p/ p + f^/2 + is the Bernoulli group of mechanical energy terms, u; = V x -u is 
the vorticity, -u x cJ is the inertial vortex force that causes turbulence, uV^v is the viscous 
force that damps it out, Fq = — V0 is the gravitational force and has been absorbed in B, 
(f) is the gravitational potential energy per unit mass in the expression = AnpG, G is 

— * 

Newton's constant, Fm is the magnetic force, and other forces have been neglected. Eq. 4 
applies in a gas or plasma when a sufficient number of particles are assembled, so that the 
particle separation Lp and the collision distance Lc are much smaller than the size L of the 
assemblage or the scale of causal connection Lh = ct, where c is the speed of light and t is 
the age of the universe. Turbulence develops whenever the inertial-vortex force of the flow is 
larger than the other terms; that is, if the Reynolds number Re = [v x uj) / {u'V'^v) , Froude 
number Fr = [v x uj)/Fg, and all other such dimensionless groups exceed critical values. 

In Landau and Lifshitz 1959 §23 an exact solution of Eq. 4 for an incompressible viscous 
fluid attributed to G. Hamel 1916 (usually termed the Jcffrcy-Hamcl flow. White 1991) gives 
multiple maxima and minima for expanding flows between plates. This solution is used 
to illustrate the relative instability of expanding flows compared to converging flows. The 
Jeffrey-Hamel converging flow solution approaches the solution for converging ideal fluid 
flow and thus might appear to be stable to the formation of turbulence at high Reynolds 
number because the turbulent intensity 5v/v decreases along a streamline as v increases. 
Converging sections are used in wind and water tunnels before test sections to decrease 
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the turbulent intensity, but the turbulent viscous dissipation rate e and turbulent velocities 
actually increase in such flows (Batchelor 1953). Can steady inviscid flows of any kind be 
stable? 

What about the stability of the expanding universe which is not incompressible but is 
a uniform expansion with rate-of-strain 7 1/t, where 7 is generally termed the "Hubble 
constant" and the expansion is termed the "Hubble flow"? Instead of decreasing along 
a streamline with 1/a: as for the incompressible diverging Jeffrey-Hamel flow, the speed 

V ^ '-fx increases with distance x. Does this mean the expanding Hubble flow is stable, 
similar to the converging JcfTrcy-Hamcl flow where the speed also increases with distance? 
Does this mean that the small CMB temperature anisotropics simply reflect the fundamental 
stability of a Hubble flow, and does not imply that large viscous or buoyancy forces must 
have been present in the plasma epoch? Are self gravitating fluids fundamentally different 
from stratified natural fluids in that the first turbulence of the Hubble flow is caused by 
gravitational forces rather than inhibited by them as in stratified flows? 

According to the further analysis and discussion in Landau and Lifshitz 1959, in §27 
titled "The onset of turbulence", steady inviscid flows are absolutely unstable. Thus, all 
flows should develop turbulence at high enough Reynolds numbers, including the diverging 
Hubble flow of the expanding universe. In their derivation, maximum amplitudes of Fourier 
modes ~ {Re ~ RccritY^'^ are expressed as functions of their departures from critical 

Reynolds numbers Rccrit and it is shown that the individual modes grow to finite values 
with increasing Reynolds number, but with an ever increasing number of modes as Re 00. 
Landau-Lifshitz admit that prediction of the mode amplitudes is mathematically difficult 
and that such stability analysis has had limited success in predicting the transition to tur- 
bulence except to confirm the 1883 Reynolds criterion, for which there is no experimental 
counterexample. As we have seen, in apparent counterexamples such as the Jeffrey-Hamel 
converging incompressible flow the increasing velocity along streamlines masks the develop- 
ing turbulence, but does not prevent it. 

The absolute instability of steady inviscid flows can be understood from the first two 
terms of Eq. 4, shown in Eq. 5. Such a flow must be irrotational to remain steady with 
dv/dt = and B constant. Otherwise the vorticity uj would produce inertial vortex forces 

V X uj that would spread the rotational region indefinitely to larger and smaller scales by 
undamped turbulent diffusion. If a variation in speed occurs along one of the steamlines, 
then accelerations 

dv/dt = -VB (5) 

develop that amplify any perturbations in v with increasing time. Increasing v requires 
increases in both B and its gradient, and decreasing v decreases both B and its gradient. 
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Prom Eq. 5, positive speed perturbations increase —VB and cause speed increases, and 
negative speed perturbations cause decreases in — V-B and cause speed decreases. Vorticity 
cJ > develops and forms turbulence, which will grow in size and kinetic energy. This 
positive feedback is independent of the continuity equation or the equation of state for the 
fluid. Finite length scale perturbations of any of the hydrophysical parameters {v,p,p) in a 
steady, inviscid, irrotational flow will cause local perturbations in the vorticity on the same 
finite scale, with resulting formation and growth of turbulent inertial vortex forces v xcu and 
thus turbulence at larger and smaller scales, drawing energy from the assumed variations 
of V along streamlines. Even the extreme case of steady flow is unstable to a vorticity 
perturbation, since the rotational region of the u) perturbation without viscous damping will 
spread its vorticity and kinetic energy, and thus turbulence, to indefinitely larger volumes 
by turbulent diffusion. 

We conclude that steady inviscid flows are absolutely unstable, confirming the 1959 
Landau- Lifshitz result and the conventional Reynolds criterion for turbulence formation. 
Viscosity is not necessary to the formation of turbulence, only its evolution. From the 
vorticity conservation equation following a fluid particle in a fluid with variable density 

DO/Dt = 80/ dt + (v ■ V)uj ^uj-i+(Vpx Vp) /p^ + uV^uj (6) 

we see variations in the density of the fluid can produce vorticity if pressure and density gra- 
dients arc not aligned, at rate (Vp x Vp)/p^j leading to unconstrained inertia! vortex forces 
vxu and thus turbulence. Vorticity is produced by vortex stretching at a rate uj ■ e, where e 
is the rate of strain tensor. Turbulence is defined as an eddylike state of fluid motion where 
the inertial vortex forces of the eddies are larger than any other forces that tend to damp the 
eddies out (Gibson 1999a). Turbulence always starts at the smallest possible scale permitted 
by viscous forces, and cascades to larger scales by a process of eddy pairing and entrainment 
by the turbulence of irrotational fluid (Gibson 1991). Fourier modal analysis fails to properly 
describe the formation of turbulence, gravitational structure formation, or small scale turbu- 
lent mixing at small Prandtl numbers. These failures result from sacrificing realistic physical 
models for mathematical convenience by considering the linear behavior of sine waves rather 
than the nonlinear behavior of finite-scale local perturbations (Gibson 1996). 

What about cosmic drag? It is sometimes argued that turbulence is prevented by the 
expansion of the universe because momentum decreases as V{t) = Vo/a{t) from general 
relativity, where Vq is an initial velocity perturbation and a{t) is the cosmic scale factor 
which monotonically increases with time t as the universe expands. Although the momentum 
and velocity of a perturbation may decrease, the proper length scale of the perturbation 
L{t) — Loa{t) will increase, so that a{t) in the Reynolds number Re{t) = VL/u will cancel. 
To first order, the Reynolds number after infiation and before mass-energy equivalence is 
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i?e 1 because V !V L ct, and u cH. Before inflation, much larger Reynolds numbers 
were possible (Gibson 2000). 

It is not true that simply because the initial perturbation of a nonlinear process is small 
that the process can be accurately described by linear theories. In particular, just because 
remnant density perturbations Sp/p ~ 10""^ from big bang quantum gravitational chaos are 
small does not mean that their evolution can be accurately described by linear methods 
once they reenter the horizon. Decreasing the size of Sp/p by a factor of 10~^ increases 
the gravitational condensation time by less than a factor of two. Cold dark matter theories 
that suggest an acoustic peak in the CMB temperature spectrum are therefore questionable. 
The gravitational response to density perturbations is always nonlinear and requires non- 
linear fluid mechanics for its description independent of the size of the perturbation or the 
predictions of the hnear, acoustic theory of Jeans 1902. 

What about energy conservation? Won't pressure support or thermal support pre- 
vent gravitational condensation at scales smaller than Lj7 Won't continued gravitational 
collapse require a loss of thermal energy to prevent pressure stabilization, and won't this 
require a spontaneous and highly efficient flow of heat from a cold object into a hot en- 
vironment? These misconceptions are all part of Jeans's 1902 legacy. Consider a volume 
of initially stagnant, constant density gas, smaller than the horizon, with mass perturba- 
tion M' suddenly placed near its center. This system is absolutely unstable to gravita- 
tional condensation or void formation, depending on whether M' is positive or negative. 
Cravitational acceleration starts immediately with radial velocity Vj. ^ —tGM'(t)/r'^, and 
mass flux Anr'^pVrit) ^ 4:7rpGM'{t)t = dM'{t)/dt independent of radius. Thus M'{t) = 
M' {0)exp{2TcpGt'^) = M'{0)exp[2n{t/TG)'^]. The density, temperature, and dynamical pres- 
sure p/p + v'^/2 remain constant during the gravitational free fall process except in the small 
space-time region of the nonacoustic density nucleus at r ^ Lj and t ^ tq = {pG)~^/'^ 
(Gibson 1999a, Gibson & Schild 1999a). Everything happens at once when t — > tq- Since it 
takes t ^ tg for information to propagate a distance Lj, no pressure support mechanism is 
possible to prevent the self gravitational collapse or void formation at nonacoustic density 
perturbations. 

In any real fluid, the Hubble flow is unstable at all scales where the Reynolds number 
exceeds a universal value Rccru ~ 100. Thus, Re ^ Sv x x/u jx^/v ~ 100 at a critical 
length scale Xcrit ~ 10(1^/7)^/^. The viscous dissipation rate e ^ 1/7^, so 



X^U ~ 10(i/3/7£V2)l/2 ^ 



(7) 



where 




(8) 
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is the 1941 Kolmogorov length scale. Turbulence always begins at scales of lOL^ and is 
inhibited at smaller scales by viscous forces. These small eddies pair, pairs of eddies pair with 
other eddy pairs, and so forth. Irrotational (and therefore nonturbulent) fluid is entrained 
into the interstices of the turbulent domain as ideal flows, is made turbulent at Kolmogorov 
scales by viscous forces, and supplies the kinetic energy of the turbulence. We now use these 
results to examine the formation of turbulence, and its inhibition, during the plasma epoch 
before 10^^ s (300,000 years). 

3. The plasma epoch 

What about the formation of turbulence in the plasma epoch? Since the Hubble flow 
is unstable to the formation of turbulence, either viscous forces or buoyancy forces, or both, 
must have been present to prevent strong turbulence. When did the flrst turbulence form? 
What was the viscosity of the plasma required to prevent turbulence? 

From COBE to WMAP, numerous experiments have been undertaken to resolve the 
small scale fluctuations of the CMB. Super-horizon contributions to the 5T variance are ap- 
proximately constant with a Sachs- Wolfe plateau of about 2 x 10~^ K for angular separations 
9 greater than about 1-2 degrees corresponding to the horizon scale Lh ~ 3 x 10^-*^ m existing 
at this plasma-gas transition time 10^^ s (Lineweaver 1999). Prom measurements at smaller 
sub-horizon scales a sonic, or doppler, peak of about 8 x 10~^ K 9 ^ 0.5 degrees and 
smaller-amplitude, smaller-scale, harmonics are attributed to undamped soTind waves in the 
plasma sloshing in CDM clump potential wells in the gravitational potential. 

This sonic peak explanation of the CMB is questionable for at least three reasons: 1. 
the postulated CDM fluid with Lsd > Lh (see Eq. 11 below) is too diffusive to condense; 
2. no sound source of any kind exists, and certainly not the non-turbulent super-powerful 
sound source that would be required to match the observations; 3. even if a super-powerful 
source of sound could be identifled, the sound would be rapidly damped by viscous forces 
because the sonic attenuation coefficient a ~ v/Vs\^ is 3> since v ~ VsLq and Lc S> A 
for all the relevant sonic wavelengths A. 

Reason 3. is why sonic fluctuations of temperature in the relatively noisy atmosphere 
of the earth rarely exceed the 1 db reference level 5T/T ^ 10^^°, Pierce and Berthelot 
1990, and why whales near Japan can be heard from California but eagles cannot. Time 
tps ^ 10^^ s (30,000 years) is indicated as the time of first structure formation since this 
is the time when the increasing horizon mass p{ct)^ just matches the observed mass of 
superclusters fa 10^^ kg (Gibson 1997b). This supercluster mass is 10~^ times the present 
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horizon mass {ct)^pcrit — 10^^ kg since the observed supervoid size is 10~^ x Lh- The 
observed globular star cluster density ~ 10~^^ kg just matches the baryonic density 
existing at t ~ 10^^ s, also indicating tps ^ 10^^ s as the time when the plasma first began 
fragmentation. Voids formed at that time should expand for a brief period as rarefaction 
waves with velocities limited by the sound speed Vs = c/S^^^, giving a structural rather than 
sonic peak in the range 0.6 > 9sp > 0.1 degrees, as observed, with a monotonic decrease of 
the 5T power spectrum reflecting fragmentation to galactic scales, possibly with acoustical 
harmonics from the rarefaction waves. Further fragmentation at smaller and smaller scales 
limits the amplitude oi 5 p/p to small values as Msv decreases toward proto-galaxy masses 
in the cooling, expanding plasma. 

To prevent turbulence at the horizon scale at decoupling requires a viscosity Ucru ~ 
cH /lOO ~ 10^^ m^ s""^, which is too large for the baryonic component by any known 
mechanism. Setting Xcrit — Lh — lOL^ in Eqs. 7 and 8 with j — 1/t gives a value of 
v = (ct/100)^7 = 9 X 10^^ m^ for our estimated tps ~ 10^^ s. This large value of u is 
only slightly larger than that required to prevent turbulence at the time of first structure. 
Once gravitational structure formation begins, buoyancy forces will inhibit turbulence. 

Densities were larger at this earlier time (30, 000 yr) so mean free paths for collisions 

Lc ~ {an)~^ were shorter, where a is the collision cross section and n is the particle density. 
The physical mechanism of viscous stress in the plasma epoch is photon collisions with the 
free electrons of the plasma (Silk & Ames 1972, Thomas 1930). The electrons then drag along 
the protons and alpha particles of the primordial plasma to maintain electrical neutrality. 
The kinematic viscosity is then 

1/ f» Lc XV — c/uTUe (9) 

where = 6.65 x 10~^^ m^ is the Thomson cross section for scattering and Ue is the number 
density of the free electrons. Substituting Ue ~ 10^° m~^ for the electron number density 
Bit t — 10^^ s (Weinberg 1972) gives u A x 10^^ m^ s~^, which is close to our estimated 
minimum v value required to inhibit turbulence. The collision distance Lc ~ 1.5 x 10^^ 
m is less than the horizon scale L^f = 3 x 10^° m, so the assumption of coUisional fluid 
dynamics in Eq. 9 is justified. The viscous dissipation rate e ~ 1/7^ 4 x 10^^ s~^ gives 
a Kolmogorov scale Lk ~ 2 x 10^° m from Eq. 8. Since lOL^ > Lh, the Hubble flow of 
plasma should be viscous and laminar or weakly turbulent. 

The baryonic density ai, t ^ 10^^ s was p ^ 2 x 10"^''' kg m~^ (Weinberg 1972). The 
strain rate at turbulence fossilization was 10~^^ s~^. Thus, from Eq. 3 

Lsv ~ (10"^^ X 4 X 10^72 X 10"^^ X 6.672 x 10"^^)^/^ ^5x 10^° m, (10) 
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approximately matching the horizon scale Lh = 3 x 10^° m. The horizon scale baryonic 
mass Mh = L\ x p = {ctYp = 5 x 10^^ kg is close to the baryonic mass of superclusters 
{Msc ~ 10^^ kg includes the non-baryonic component), so these are suggested as the first 
structures of the universe, formed by fragmentation when the viscous Schwarz scale first 
matches the horizon scale, Gibson 1996, 1997ab, 1999ab. 

Proto-superclusters formed by fragmentation rather than condensation because void 
formation is augmented by the expansion of the universe but condensation is inhibited. Thus 
proto-superclustcr-voids expand in the plasma epoch while the proto-superclusters between 
these voids also grow, but more slowly, by internal fragmentation, preserving the density of 
the fragments. Further fragmentation at Lsv scales down to protogalaxy masses with little 
change in the baryonic density due to fossil density turbulence formation is proposed by 
HGT (Gibson 1996). Turbulence formation is inhibited at every stage of the plasma epoch 
by a combination of viscous and buoyancy forces, and there is no energy source for sound 
other than the gravitational void formation. Temperature fluctuations observed in the CMB 
are proposed as fossils of big bang turbulence and fossils of the first structure formation, 
Gibson 2000. 

In contrast. Silk 1989 Fig. 10.1 traces the evolution of an adiabatic galaxy mass pressure 
fluctuation as it drops below the Jeans mass at a redshift of z ^ 10^ and oscillates as an 
undamped sound wave in the necessarily inviscid plasma epoch with 10"^ density contrast until 
decoupling at z ^ 10^. It seems unlikely that any such loud sounds (> 100 db) could start at 
that time, less than a week after the big bang. If somehow they were started they would be 
rapidly damped, within another week, by the large photon viscosity, not to mention damping 
by the expansion of the universe (cosmic drag). Sonic pressure fluctuations p ~ Poexp[—ax], 
where Po is the initial pressure, x is the direction of propagation. The sonic attenuation 
coefficient a ^ uuj'^/Vg — u/VsX^, where the frequency u; — V5/A, A is the wavelength, and 
Vs is the sound speed (Pierce and Berthelot 1990). Thus, p/po ~ exp[—{iy/VsX'^)x] <^ 1 for 
distance x > X if u > VsX, and this will be true since i' increases with time as the universe 
density decreases and X < ct is limited in size by the time t when the sound wave was 
created. 

What about the non-baryonic dark matter (NB) required to make up the critical density 
of a flat universe? Its cross section ac for collisions with ordinary matter must be very small 
or it would have been detected based on the expression ac = rrip {GM/ry/'^/pDNB, where 
rrip is the particle mass and -Datb is the diffusivity inferred from outer-halo dimensions r of 
galaxies or clusters of mass M (Gibson 2000). Thus such material must have large mean free 
paths for collisions and large diffusivities D^b = L^b x vnb compared to for baryonic 
matter since Ljvb ^ Lb and vnb ~ vb- From measurements of the mass profile of Abell 
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1689 by Tyson & Fischer 1995, Gibson 1999b estimates the non-baryonic dark matter of the 
dense galaxy cluster is /^atb ~ 10^^ m^ s~^ by setting the radius of curvature of the profile 
to LsD- Since neutrinos are now known to have mass, an obvious non-baryonic candidate is 
neutrinos, which have densities comparable to the density of photons and very small cross 
sections for collisions since they interact with baryonic matter mostly through the weak force. 
Large numbers of neutrinos were formed in nucleosynthesis and their unknown number of 
flavors and abihties to convert between flavors leaves their total mass a mystery. Assuming a 
neutrino coUision cross section of (7„ ^ 10~^° m^ and number density n„ ^ 10^° m~^ gives a 
mean free path L„ of 10^° m, so coUisional dynamics apply. Cross sections for light (fs 10~^^ 
kg) particles like neutrinos give such small a values, but 10~^^ kg particles like neutralinos 
give a ^ 10~^^ m^, much larger than ^ 10~^^ m^ theoretical values or the < 10~^^ m^ values 
excluded by experiments (Gibson 2000). 

In the case of strongly diffusive matter in weakly turbulent flows, gravitational conden- 
sation is limited by a match between the diffusion velocity of an isodcnsity surface Vd ~ D/L 
and the gravitational free fall velocity Vq ~ L/tq, giving the diffusive Schwarz scale 



where the diffusivity Dn ~ L„ x c ~ 3 x 10^^ m^ s~^. This gives Lsd ~ 10^^ m during the 
plasma epoch, much larger than any of the structures formed and larger than the horizon for 
part of the epoch. Any such nonbaryonic material would diffuse away from the protogalaxies 
and proto-superclusters as they fragment, to fill the voids between. Non-baryonic materials 
fragment as the last stage of gravitational structure formation to form protosuperhalos when 
the baryonic protosuperclusters separate by scales larger than Lsd- This is contrary to cold 
dark matter models that require CDM condensation as the first rather than last stage of 
structure formation, producing, rather than being produced by, the baryonic structure. 

The necessary condition for the diffusive Schwarz scale Lsd of Eq. 11 to determine the 
minimum scale of gravitational condensation is 



for viscous flows. Since ^tq > 1 and D ^ v for baryonic matter, the scale Lsd only applies 
to nonbaryonic matter. Substituting D„ 3 x 10^^ m^ s~^ and p 10~^^ kg m~^ for 
the density of a galaxy cluster gives Lg^i 3 x 10^^ m (Mpc) as the scale for gravitational 
fragmentation of the non-baryonic dark-matter halo of a small galaxy cluster with total mass 



Lsd ^ {D^lpG) 



1/4 




D > v^tq 



(12) 



4 X 10^^ kg. 



Thus a proper description of structure formation in the primordial self-gravitational 
fluids of the early universe requires more than the linearized Euler equation with gravity 
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and the density equation without diffusion or gravity, as assumed by (Jeans 1902). All the 
forces in the momentum Eq. 4 are needed except [Fm + •••)• The appropriate non-acoustic 
density conservation equation near density maxima and minima is 

dp/dt + v-Vp = DeffV^p, (13) 

where -De// = D — L'^/tq and D is the molecular diffusivity of the density; that is, on scales 
L < Lj so that the pressure adjusts rapidly, and on scales L > Lsxmax (the maximum 
Schwarz scale) where gravity dominates fluid forces and molecular diffusion (Gibson 1999b). 
The density p depends on temperature and species concentration variations and their dif- 
fusivities, and not simply the pressure as assumed by Jeans. The problem is similar to the 
turbulent mixing problem (Gibson 1968) except for the remarkable fact that for clouds of 
fluid with sizes Lj > L > Lsxmaxi gravitational diffusivity takes over and the effective diffu- 
sivity Df.ff becomes negative. Thus, rather than reaching a local equilibrium between local 
straining and diffusion at the Batchelor length scale Lb = {D/'jy^'^ near density extrema as 
in turbulent mixing theory with a monotonic decrease toward ambient values, densities in 
the self-gravitational fluids of astrophysics increase to large values or decrease toward zero 
at these points due to gravitational instability (Gibson and Schild 1999a). 



4. The gas epoch 

From standard cosmology and the CMB observations, the initial conditions of the gas 
epoch are precisely defined. Little or no turbulence was present, as discussed previously, so 
the rate of strain of the fluid was larger than 7 ^ 1/t ^ 10^^^ s^-"^ existing at that time 
and smaller than the fossil vorticity turbulence value in the structures '-fps ~ 10^^^ s^^. 
The density of the protogalaxies cannot have been much different from the fossilized initial 
fragmentation density pps ~ 10~^^ kg m~^ since there was insufficient time for collapse. The 
temperature at decouphng was ^ 3 000 K. The composition was 75% H and 25% He by 
mass. Therefore the kinematic viscosity of the primordial gas was about 3 x 10^^ m^ s~^, from 
// = p X 1/ in standard gas tables with a weighted average /x(T'o), with gas constant R about 
3 612 m^ s~^ K~^; (Gibson 1999b). Viscous dissipation rates were only e ~ ?a 3 x 10~^"^ 
m^ s~^ so the Kolmogorov scale Lx ~ 5 x 10^^ m from Eq. 10 was a factor of 5 x 10^ smaller 
than the horizon scale. 

Fragmentation of the neutral gas protogalaxies occurred simultaneously at both the 
Jeans scale Lj of Eq. 1 and the viscous Schwarz scale Lsv of Eq. 3, where for the primordial 
gas conditions Lj ^ lO^L^y ^ Lsv- The physical mechanism of this Jeans scale fragmenta- 
tion is not the mechanism proposed by Jeans 1902. Temperatures in growing voids at scales 
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smaller than Lj ^ {RT / pGY^'^ adjust by particle diffusion to remain constant at T = p/pR 
as the gravity driven rarefaction waves of void formation propagate, where the term "void" 
indicates a density deficiency rather than p = 0. As the density decreases the pressure de- 
creases. Particle speeds and temperatures are constant as long as the particle diffusion time 
Tp = L/ [RTy/"^ is less than the gravitational free-fall time r^; that is, for scales L < Lj. For 
scales L larger than Lj the diffusion time Tp is larger than tg, causing temperatures in these 
large voids to decrease as the voids grow because particle diffusion cannot maintain constant 
temperature and acoustical equilibrium. When this happens, radiation heat transfer from 
the warmer surroundings increases the temperature, and thus also the pressure, within the 
voids, and the increased pressure accelerates the void formation, isolating blobs of gas at 
some multiple of the Jeans scale to form PGCs. 

Substituting the T and p values of the primordial gas gives Lj fa 5 x 10^^ m, and 
Mj = Ljp ^ 10^5-36 Substituting p, 7 and ly values in Lsv = {i^l/pGY^'^ gives 
Lsv ~ 10^'' m and Msv = LlvP ~ lO^^'-^s kg, a factor of pa 10^^ smaller than Mj. The 
Jeans scale objects are called "Proto- Globular- Clusters" (PGCs) and the Lsv scale objects 
are called "Primordial Fog Particles" (PFPs). From the observational evidence it appears 
that many if not most PGCs have not yet dispersed and most of their PFPs have not yet 
accreted to form stars, so that both persist as the dominant component of galactic baryonic 
dark matter (Gibson 1996). The calculated masses of PGCs and PFPs depend on univer- 
sal proportionality constants of order one that will emerge from observations. Observations 
of globular star clusters indicate a mass 10^~^Mq matching our calculated PGC value of 
-|^q35-36 jj.g densities close to the fossilized initial fragmentation density pps ~ 10~^^ 
kg m~^. The calculated PFP mass 10^^~^^ kg matches observations of ~ 1Q~^Mq "rogue 
planets" by Schild 1996 as the dominant component of the Icnsing galaxy in a lensed quasar 
system. Evidence supporting HGT has recently been summarized (Gibson & Schild 2003), 
and includes the appearance of PFP candidates brought out of cold storage by evaporation 
near hot objects such as white dwarfs in planetary nebula. Figure 2 shows PFP candidates 
in the Hehx planetary nebula which is the one closest to earth, photographed by the Hubble 
Space telescope. Thousands of cometary globules appear with mass values, densities, and 
separation distances as predicted by HGT. 



5. Conclusions 

We conclude that the small amplitude 6T /T a; 10~^ of measured temperature fluctua- 
tions in the cosmic background radiation is evidence of strong turbulence damping by both a 
photon viscosity v ^ Ax 10^^ m^ s~^ and buoyancy forces of viscous-gravitational structure 
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formation beginning approximately 30, 000 years after the Big Bang. 

The hypothesis is rejected that the small CMB fluctuations reflect hydrodynamic sta- 
bility of the Hubble flow. This would require a critical Hubble flow Reynolds number of 
Rccth ~ 10^, contrary to universal similarity hypotheses of Kolmogorov 1941 for turbulence 
and strong experimental evidence that the universal critical Reynolds number of transition 
is Rccr < 25. Steady inviscid flows arc absolutely unstable to the formation of turbulence, 
as shown in §2 and as derived by Landau & Lifshitz 1959. Buoyancy forces from gravita- 
tional structure formation in the plasma epoch are therefore required to explain the lack of 
turbulence at the time of plasma to gas transition. 

The hypothesis is rejected that the CMB spectral peak with L ^ OAL^ is a dopplcr 
or sound horizon with L ^ Vst because no source of sound exists to produce the observed 
6T/T ^ 8 X 10^^ peak value, and because strong viscous damping in the plasma epoch 
would rapidly flatten any such sonic peaks. Persistent sonic oscillations of baryonic matter 
sloshing in CDM potential wells as a sound source is rejected because no CDM potential 
wells are possible in the plasma epoch, because viscous damping would occur, and because 
recent strong observational evidence excludes CDMHCC scenarios, Gibson & Schild 2003. 
Instead, the observed spectral peak at scales L ct is interpreted as evidence of the flrst 
hydro-gravitational structure formation. Secondary acoustic peaks observed may reflect rar- 
efaction wave oscillations of hydro-gravitationally driven proto-supercluster void formation 
near sonic velocities. Hydro-gravitational theory suggests the first structures to form were 
proto-supcrclustcr- voids at the viscous Schwarz scale Lsv, when Lsv > Lh first matched 
the increasing horizon scale Lh- Rapid expansion of the universe during the plasma epoch 
prevented gravitational condensation but enhanced void formation, §3. 

As shown in §4, fragmentation of the primordial gas occurred simultaneously at Lj 
and Lsv scales to form proto-globular-clustcrs (PGCs) and primordial- fog-particles (PFPs). 
Estimated PGC masses match the observed globular star cluster masses of lO^M© and es- 
timated PFP masses match the observed "rogue planet" dark matter masses of 1O~^M0 
in lensed quasars (Schild 1996). Most PGCs and their PFPs are observed to persist as 
dark clumps of frozen planetoids, forming the dominant component of ^ lOOkpc — 3 x 10^^ 
m galactic dark- matter inner-halos (Gibson & Schild 2003), with non-baryonic dark matter 
fragmenting to form outer galactic dark-matter halos at Lsd ~ Mpc — 3 x 10^^ m scales. 

All evidence suggests the early universe was an extremely gentle place, with practically 
no turbulence or sound anywhere after the big bang and prior to the formation of stars. 

Biioyancy and large photon viscosities damped sound and turbulence in the plasma epoch. 
Gravitational condensation formed PFPs and prevented turbulence in the early stages of the 
gas epoch. As the universe continued to expand and cool, some of these small-planetary-mass 
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objects experienced an accretional cascade to larger mass-scales to form the first very small 
stars. This cascade was a gentle process to produce the remarkably spherical distributions of 
these long-lived, tightly-packed stars in globular star clusters and the more numerous dark 
or dim PGC-PFP baryonic-dark-mattcr structures with the same p ^ 10^^^ kgm^'^ density 
and the PFP mass ~ 10^^ kg preserved as fossils of the weak turbulence and large density 
at the time of first structure (Gibson 2000). 

The author is grateful for many constructive questions and comments from Rudolph 
Schild. This paper is dedicated to the memory of George Keith Batchelor, 8 March 1920 — 
30 March 2000. 
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V 

Fig. 1. — Turbulent Hubble flow wakes of proto-galaxies (labeled PG) bound by gravity (dark 
double arrows) in a proto-cluster, soon after the plasma to gas transition at time t ~ 10^^ s 
(300,000 years) after the big bang. From HGT (Gibson 1996) the PG size is ~ 10^° meters, 
with primordial H-He mass 10^^ kg (a; IO^^Mq) and density p fa 10~^^ kg m~^. The 
Hubble flow velocity vh = r x^h ^ 10^ m at the radius r of the PGs. Hubble flow drag 
forces will separate the protogalaxies, as it separated the protoclusters and protosuperclusters 
also formed in the plasma epoch by gravitational fragmentation. Nonbaryonic dark matter 
moves freely through the galaxies by diffusion and the Hubble flow to flll the voids. The 
horizon scale Lh where vh — cis?) x 10^^ m (100 kpc), 10 times the size of the outer sphere 
shown. 



Fig. 2.— PFP-like objects (Gibson & Schild 2003) observed by HST in the Helix Planetary 
Nebula, 4.5 x 10^^ m from Earth (O'Dell & Handron 1996). Evaporation of the frozen H-He 
objects produces 10^^ kg photo-ionized cocoons, with H-wakes pointing away from the hot 

50,000i^') White Dwarf. The indicated density of the PN halo is puaio ~ MpFpL'l^ 
10"^'' kg m~^, where the PFP separation distance Lsep. ~ 10^^ m. This matches the baryonic 
density at the time of first structure 10^^ s (30,000 years) as a fossil of this time of first 
structure formation as predicted by HGT. 
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Table 1. Length scales of self-gravitational structure formation 



Length scale name 


Symbol 


Definition" 


Physical significance 


Jeans Acoustic 


Lj 




ideal gas pressure equilibration 


Chandrasekhar Turbulent 


LjCT 


[{p + PT)/p'G]y' 


turbulence balances gravitation 


Schwarz Diffusive 


LsD 




Vd balances Vq 


Schwarz Viscous 


Lsv 




viscous force balances gravitational force 


Schwarz Turbulent 


LsT 




turbulence balances gravitation 


Kolmogorov Viscous 


Lk 




turbulence force balances viscous force 


Batchelor Diffusive 


Lb 




diffusion balances strain rate 


CoUision 


Lc 




distance between particle coUisions 


Horizon, Hubble 


Lh 


ct 


maximum scale of causal connection 



'^Vs is sound speed, p is density, G is Newton's constant, p is pressure, pt = plSv)"^, v is velocity, 
D is the diffusivity, Vd = D/L is the diffusive velocity at scale L, Vq = L[pG]^/^ is the gravitational 
velocity, 7 is the strain rate, 1/ is the kinematic viscosity, e is the viscous dissipation rate, m is the 
particle mass, a is the collision cross section, light speed c, age of universe t. 

•^Magnetic and other forces (besides viscous and turbulence) are negligible for the epoch of primordial 
self-gravitational structure formation considered here (Gibson 1996). 



